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ABSTRACT 

We present results from quantitative modeling and spectral analysis of the high mass X-ray binary 
system Vela X-l obtained with the Chandra High Energy Transmission Grating Spectrometer. The 
observations cover three orbital phase ranges within a single binary orbit. The spectra exhibit emission 
lines from H-like and He-like ions driven by photoionization, as well as fluorescent emission lines 
from several elements in lower charge states. The properties of these X-ray lines are measured with 
the highest accuracy to date. In order to interpret and make full use of the high-quality data, we 
have developed a simulator, which calculates the ionization and thermal structure of a stellar wind 
photoionized by an X-ray source, and performs Monte Carlo simulations of X-ray photons propagating 
through the wind. The emergent spectra are then computed as a function of the viewing angle 
accurately accounting for photon transport in three dimensions including dynamics. From comparisons 
of the observed spectra with results from the simulator, we are able to find the ionization structure 
and the geometrical distribution of material in the stellar wind of Vela X-l that can reproduce the 
observed spectral line intensities and continuum shapes at different orbital phases remarkably well. 
We find that the stellar wind profile can be represented by a CAK-model with a star mass loss rate 
of (1.5-2.0) xlO _6 M0yr -1 , assuming a terminal velocity of 1100 km s _1 . It is found that a large 
fraction of X-ray emission lines from highly ionized ions are formed in the region between the neutron 
star and the companion star. We also find that the fluorescent X-ray lines must be produced in at 
least three distinct regions - (1) the extended stellar wind, (2) reflection off the stellar photosphere, 
and (3) in a distribution of dense material partially covering and possibly trailing the neutron star, 
which may be associated with an accretion wake. Finally, from detailed analysis of the emission line 
profiles, we demonstrate that the stellar wind dynamics is affected by X-ray photoionization. 
Subject headings: X-rays: binaries - - X-rays: individual (Vela X-l) - - stars: neutron - 
starts: winds, outflows 



1. INTRODUCTION 

In a high mass X-ray binary (HMXB), a compact ob- 
ject, either a neutron star or a black hole, sweeps up 
material as it orbits through the stellar wind of a mas- 
sive O- or B-type companion star. When the material 
is accreted onto the compact object, a fraction of its 
gravitational potential energy is converted into X-rays, 
which in turn ionizes and heats the surrounding stellar 
wind. X-ray photons from the compact object are re- 
processed by the ionized material, resulting in discrete 
spectral features, which carry a wealth of information 
about the geometry and physical state of the material in 
the system. The compact object can, therefore, be used 
as a radiation source to probe the structure of the stellar 
wind to derive the physical parameters that characterize 
its nature. 

Vela X-l is the archetype HMXB consisting of a neu- 
tron star and a massive companion star, and as such it 
has been the most well-studied object in this class since 
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its discovery. It is an eclipsing HMXB system containing 
a pul sar with a pulse period of 283 s iMcClint ock et al 
1976) and an orbital period of 8.964 days l)Forman et al 



The optical co mpanion star, HP 77581 
is a B0.5Ib supergiant llBrucato fc Kristianl Il972t 
iHiltner. Werner fc Osmerl Il972|) . which drives a stel- 
lar wind with a mass- loss rate in the range (1 - 
7) x 10~ 6 Mr* vr" 1 !lHutchingsll97"6HDupree et alJll98Ct 
i Kallman fc Whitell982USadakane et alJl985HSato et alJ 
Il986|) . The terminal velocity of the stellar wind was de- 
termined to be 1100 k m s" 1 from the P -Cygni profile 
of UV resonance lines l)Prinia et alJ 199Q 1 ). Its average 
intrinsic X-ray luminosity is ~ 10 36 erg s _1 , which is 
consistent with accretion of the stellar wind by the grav- 
itational field of a neutron star. 

The first attempt to model the global X-r ay emission 
line sp ectrum of Vela X-l was presented by iSako et all 
(1999) using the CCD spectrum obtained with ASCA. 
They characterized the wind u sing a standard velocity 
profile of OB stars proposed bv lCastor. Abbott fc Klehl 
(1975) (CAK model). By adjusting the mass loss rate, 
they found a statistically acceptable fit to the ASCA spec- 
trum. However, since the energy resolution of ASCA is 
not sufficiently high to resolve any of the emission lines 
predicted by the ionization model, the intensity measure- 
ment of each line is limited due to blending. In addition, 
since the expected X-ray line shifts and widths are much 
smaller than the width of the instrument response, they 
were not able to extract any information about the dy- 
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namics of the X-ray emitting gas. 

The spectral resolution (E/AE ~ 100-1000) of the 
transmission grating spectrometers on board Chandra 
have provided a detailed view of the emission line 
spectrum, including accurate measurements of line in- 
tensities an d velocity profiles i n a number of HMXBs 
(Vela X-l: lSchulz"etaTI 1200% iGoldstein etafl 12001 
GX3 01-2: IWatanabe et al.l2003t Cyg X- 3: lPaerels et al.1 
2000; 4U 1700-37: IBoroson~ ai] l20f)l . ISchulz et al.1 
( 2002) present the first high-resolution X-ray spectrum of 
Vela X-l observed during eclipse. They detect a narrow 
RRC from Ne X with an electron temperature of ~ 10 eV, 
which implies that the stellar wind is in deed driven 
by photoionization. iGoldstein et alJ l)2004|) present the 
spectra of Vela X-l observed during the three orbital 
phase ranges - eclipse, <j> = 0.25, and <f> = 0.50. They 
report the first detection of multiple emission lines at or- 
bital phase = 0.5 and analyzed the spectral changes 
between the three orbital phases. 

In this paper, we present a quantitative analysis 
of the X-ray s p ectra of Vela X-l first presented by 
IGoldste m "etall l|2"004|) using the data obtained with 
Chandra High Energy Transmission Grating Spectrom- 
eter (HETGS) during three orbital phase ranges. We 
construct a simulator to compute the physical condi- 
tions of the surrounding gas for a given binary config- 
uration, stellar wind parameters, and X-ray luminos- 
ity. Continuum photons that originate from the com- 
pact object interact with the atoms to produce emission 
lines and continua, which are propagated through the 
wind accounting for all relevant radiative transfer effects 
including line scattering/absorption, continuum absorp- 
tion, and Compton scattering. The emerging spectra 
are calculated and compared with the observed orbital- 
phase-resolved spectra, which provide constraints on the 
global physical properties of the ionized stellar wind. 

This paper is organized as follows. In § 2, we describe 
the Chandra observations and summarize the data re- 
duction procedures. The reduced spectra are fitted using 
simple empirical models to extract the general properties 
of the stellar wind. In § 3, we provide a full description 
of the simulator, which is used to model the binary sys- 
tem and to calculate the resulting spectra as a function 
of the orbital phase. In § 4, we describe a self-consistent 
picture of the Vela X-l system derived from the simula- 
tor in comparison to the data. We finally summarize our 
results in § 5. 

2. OBSERVATIONS AND RESULTS 

Vela X-l was observed with the Chandra 
HETGS/ACIS-S during three different orbital phase 
ranges within a single binary orbit in February 2001 
- (1) <f> = 0.237-0.278, (2) = 0.481-0.522, and (3) 
4> = 0.980-0.093, hereafter referred to as phase 0.25, 
phase 0.50, and eclipse, respectively. The observation 
dates and exposure times are summarized in Table ^ 
A schematic drawing of the binary orbit and the phase 
ranges are shown in Fig. ^ 

All of the data are processed using CIAO v2.3 6 , and 
spectral analyses are performed using XSPEC 7 . Since the 
zeroth order image was severely piled-up during phase 

6 http://cxc.harvard.edu/ciao/ 

7 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 



0.25 and phase 0.50, the locations of the zeroth order 
image were determined by finding the intersection of the 
dispersed events as well as the streak events, which are 
produced during CCD frame transfer. This gives us a ze- 
roth order location that is accurate to within 0.1 pixels, 
which is much smaller than the instrument line response. 
The positive and negative spectral orders are inspected 
separately to check that systematic offsets do not exist. 
We then apply spatial filters for both the Medium En- 
ergy Grating (MEG) and the High Energy Grating (HEG) 
data, and subsequently use an order sorting mask by us- 
ing the energy information obtained with ACIS-S. Only 
the first order events are used in our spectral analyses. 
The background events arc estimated from the adjacent 
region to the dispersed event region and are found to 
contribute at most 5 % for the eclipse data and 3% for 
the phase 0.25 and phase 0.50 data. The background is 
subtracted in the spectral analyses. 

2.1. Continuum Emission 

Figure shows the HEG spectra for the three orbital 
phases. X-rays emitted from the neutron star are di- 
rectly observed during phase 0.25 and phase 0.50 as con- 
tinuum spectra. As expected from the geometry of the 
binary system, the spectrum taken in eclipse is domi- 
nated by line emission and scattered components. In 
order to parameterize the properties of the continuum, 
we fit the spectra with a photo-absorbed power-law func- 
tion. Since the spectrum of phase 0.25 is affected less by 
absorption, we leave both the hydrogen column density 
and the photon index as free parameters in the fit. On 
the other hand, for phase 0.5, we fix the photon index 
to the value obtained from phase 0.25 and calculate the 
absorption. We assume that the chemical abundance of 
metals in th e absorbing ma terial is 0.75 times the cosmic 
abundance ljFeldmad ll9921 . which is known to be a rep - 
resentative value for typical OB-stars l)Bord et al.lll976|) . 
The derived parameters from the spectral fits are listed 
in Table The best-fit models are superimposed on the 
spectra in Fig. El The absorption-corrected luminosity is 
determined to be identical for observations in phase 0.25 
and in phase 0.50 and corresponds to 1.6 x 10 36 erg s _1 
in the 0.5-10 keV ran ge, assuming a distance of 1.9 kpc 
(jSadakane et al.lll985|) . 

2.2. Line Emission 

A large number of emission lines are clearly detected in 
the spectra of phase 0.50 and eclipse as shown in Fig. OH 
and Fig. Emission lines from H-like and He-like S, Si, 
Mg and Ne can be seen, in addition to fluorescent lines 
from Fe, Ca, S, and Si ions in lower charge states. We 
note that emission lines from the same ions are detected 
in the spectra of both phase 0.50 and eclipse. 

In Fig. EI we show the three spectra in a narrow energy 
range of E = 1.70—2.05 keV, which contains the Ka tran- 
sitions of Si. An intense Lya line from H-like Si and fully- 
resolved He-like triplet lines are clearly seen in the phase 
0.50 and eclipse spectra. At slightly lower energy, a spec- 
tacular complex of emission lines are also detected, which 
we identify as Ka fluorescent lines from Si in lower charge 
states. In the phase 0.5 spectrum, the near-neutral lines 
(Si II - Si VI) are blended at E ~ 1.74 keV, and a series 
of lines above this energy are resolved into those from 
Si VII - Si XI. This implies that the stellar wind in this 
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system contains an extremely wide range in ionization 
parameter. The eclipse spectrum also shows a similar 
structure, but with much lower intensities. 

The centroid of energies, the widths and the intensities 
of each line are determined by representing them with 
single gaussian components. In the spectral fittings, we 
use the Poisson likelihood statistics, instead of the % 2 
statistics, because the numbers of photons in some of 
the bins in the spectra are small. The derived parame- 
ters for phase 0.50 and eclipse are listed in Table |3 and 
Table respectively. We also show list in Table [5] the 
line intensity ratios between phase 0.50 and eclipse for 
the H- and He-like lines that are detected with a statis- 
tical significance of more than 5 a. The H-like lines are 
brighter during phase 0.5 by a factor of 8-10 compared 
to during eclipse, and similarly the He-like lines are 4-7 
times stronger. 

The resolving power of the HETGS allows us to mea- 
sure Doppler shifts to an accuracy of ~ 100 km s _1 . 
Figure HH compares the line profiles of Si Lya and Mg 
Lya between phase 0.50 and eclipse. The offset in the 
line center energies are clearly seen. We measure the ve- 
locity shifts of all of the Ka lines of H-like and He-like 
ions, which are plotted in Fig. [7| Though some fluctua- 
tions are seen, there is a general trend that the lines are 
blue-shifted during phase 0.50, and they are red-shifted 
during eclipse as expected from an outflowing wind ir- 
radiated by an X-ray source. The velocity shifts have 
values in the range in Av ~ 300-600 km s" 1 (Table 
and widths of a < 300 km s _1 . 

Radiative recombination continuum (RRC) is detected 
clearly from H-like Ne. We fitted the RRC spectra us- 
ing the "redge" model in XSPEC. The electron tem- 
peratures are determined to be kT e = eV and 

kT e = 6.6+^8 eV during phase 0.50 and eclipse, respec- 
tively. 

Iron Ka fluorescent lines are detected in all three or- 
bital phases. The profiles of these lines are shown in 
Fig. |S1 The parameters derived from the spectral fits 
with single gaussian models are listed in Table |SJ The 
equivalent width of the iron Ka line is measured to be 
116 eV and 51 eV for phases 0.50 and 0.25, respectively. 
During eclipse, a high equivalent width of 844 eV is ob- 
served. As shown in Fig. [SJ there is an indication of a 
Compton shoulder in the iron Ka line at phase 0.50. 

2.3. Summary of Results and their Implications 

Some general global properties of the circumsource 
medium can be deduced from t he r esults of t he simple 
fitting procedures described in § 12. II and § 12.21 

The continuum from the neutron star is highly ab- 
sorbed during phase 0.5 whereas the absorption column 
density is reduced during phase 0.25. This implies that 
there is an absorber located behind the neutron star as 
viewed from the companion star. The average contin- 
uum intensities are nearly identical during the two phase 
ranges. The presence of strong soft X-ray lines during 
phase 0.5, however, suggests that this absorber is local- 
ized and does not cover the X-ray emission line regions, 
which most likely extend beyond the size of the absorber. 

The emission lines during phase 0.50 is approximately 
an order of magnitude brighter than those during eclipse. 
This indicates that a significant fraction of the X-ray line 
emission is produced in a region between the neutron 



star and the companion star, which is occulted during 
eclipse. The velocity shifts and widths provide addi- 
tional constraints on the properties of the emission re- 
gion. The measured Doppler shifts of < 500 km s _1 
are much smaller than th e terminal velocity of Vela 
X-l (v„„ = 1100 km s" 1 : iPrinia et alJ IT59TO . For a 
iCastor. Abbott fc Kleinl l)1975|) stellar wind velocity pro- 
file of isolated OB stars, this velocity is achieved at only 
0.25 stellar radii from the photosphere. The Doppler 
broadening of emission lines is also small (< 300 km s _1 ). 
These observational results shows that the emission site 
of lines from highly ionized ions is distributed on a plane 
between the neutron star and the companion star, with 
most of the emission coming from near the line of centers. 

Finally, the detection of a narrow radiative recombi- 
nation continuum from H-like Ne during phase 0.5 and 
in eclipse with electron temperatures of 6-8 eV provides 
direct evidence that the emission lines from highly ion- 
ized gas in Vela X-l are driven through photoionization. 
Other possible sources of mechanical heating are, there- 
fore, not important in characterizing the ionization struc- 
ture of at least the highly-ionized regions of the stellar 
wind. 

3. OUTLINE OF THE SIMULATOR 

3.1. Modeling of Photoionized Plasma 

The X-ray emission lines detected from Vela X-l can 
be interpreted as due to emission from a gas photoion- 
ized by the continuum radiation from the neutron star. 
When line photons are created through recombination 
cascades following photoionization, they interact with 
the surrounding gas as they propagate out of the binary 
system, and this is controlled by the ionization struc- 
ture and the density distribution of the gas in the stellar 
wind. Therefore, by investigating the properties of the 
emission lines, we can infer the physical conditions, ge- 
ometrical distribution, and the velocity structure of the 
wind. 

In order to compute the ionization structure and to 
reproduce the wealth of information provided by the 
high-resolution HETGS spectra, it is necessary to model 
the emission using realistic physical models of the stel- 
lar wind in a binary system. We have developed three- 
dimensional Monte Carlo radiative transfer code to sim- 
ulate the interactions of X-rays with the particles in the 
stellar wind. 

Our simulation code consists of two main components: 

(1) Calculation of the ionization structure. 

(2) Monte Carlo simulation for radiative transfer. 

In part (1), we construct a map of the ion abundance 
and temperature in the stellar wind. In part (2), we use 
the results from part (1) and propagate X-ray photons 
from the neutron star out of the binary system. The sim- 
ulation is performed on a three-dimensional grid of cells, 
which is capable of treating complex geometries. Radia- 
tive transfer effects that are important for highly ionized 
media are included. The medium is allowed to have ar- 
bitrary velocities, and Lorentz transformations are per- 
formed between the cells, which is necessary to calculate 
the correct output line profiles. Our detailed procedures 
are described in the following subsections. 
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3.2. Calculation of the Distribution of Ionization Degree 

The ionization balance and electron temperature are 
determined primarily by the particle density and the lo- 
cal flux and spectral shape of the ionizing X-rays. We 
compute the three-dimensional ionization structure by 
first dividing the binary system into square cells of equal 
size. Since the system is assumed to be cylindrically 
symmetric along the line of centers, we compute the ion 
fractions and temperatures on a two-dimensional grid as 
shown in Fig. Each cell is 4 x 10 11 cm along the side 
and the system consists of 61 x 31 cells. 

The density structure around the companion star is 
specified from the velocity profile determined from the 
UV observation and an assumed mass-loss rate, which 
is treated as a free parameter. This assigns a particle 
density to each of the cells. The intrinsic luminosity and 
spectrum shape of the X-ray source as determined from 
the observed spectrum are then placed in the cell where 
the neutron star is located. The ion fractions and the 
electron temperatures of the cells immediately surround- 
ing the neutron star are computed using this initial con- 
tinuum shape as the input ionizing radiation field. The 
computation moves outwards to the next layer of cells, 
this time, using a different continuum shape, which ac- 
counts for absorption by the inner cells. This procedure 
continues until the outermost cells are reached. We have 
use the XSTAR 8 program to calculate the ion fractions 
H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, and Fe, and 
the electron temperature in each cell. A flowchart of the 
procedures is shown in Fig. |5| 

This procedure ignores the diffuse radiation field as the 
source of additional ionization, which is a small contrib- 
utor to the total photoionization rate at least for the ions 
observed in our spectra. A more accurate treatment of 
this effect might involve an iterative scheme, which is not 
necessary for our present analysis, but may be important 
for extremely optically thick media. 

3.3. Monte Carlo Calculation of the X-ray Emission 
from Photoionization Equilibrium State 

The map of the ion fractions calculated by the method 
described in the previous subsection assures equilibrium 
between photoionization and recombination. We then 
use this ionization structure to propagate again the con- 
tinuum photons originating from the neutron star to 
compute the emitted X-ray spectra as a function of the 
orbital phase. This is done by accounting for radiative 
transfer effects in detail and using Monte Carlo methods 
as described below. 

The cells are set up somewhat differently from those 
used in the calculation of the ionization structure. The 
cells sizes are chosen so that each of them are not op- 
tically thick to the strongest resonance lines. Smaller 
cells are required near the companion star where both 
the density and density gradient are high, and we choose 
cubic cell sizes of 4 x 10 11 cm on each side. The cells 
are increased in size outwards and the outermost cells of 
the simulation space are larger by a factor of 7.5 com- 
pared to the inner ones as shown in Fig. ^| We start 
the Monte Carlo simulation with a photon at the position 



of the central neutron star. The photon energy is drawn 
from the observed power-law spectral energy distribution 
and is initially given a random direction. This photon is 
then traced through the cells and interacts with the ions 
in the stellar wind, which then generates a secondary 
photon through radiative recombination, radiative tran- 
sitions, and fluorescent emission. All of the photons in- 
cluding those from the central continuum and subsequent 
line photons produced by interactions are propagated un- 
til they either completely escape the simulation space or 
are destroyed by one of the physical processes described 
below. The emergent photons are then selected depend- 
ing on the viewing angle of the system, and are finally 
histogrammed to produce a spectrum. 

When a photon interacts with a particle, one of three 
physical processes is allowed to occur: (1) photoioniza- 
tion, (2) photoexcitation, and (3) Compton scattering. 
For a given photon energy, the absorption and scattering 
cross sections and hence the optical depths are calculated 
for each ion within a given cell. Interaction probabil- 
ities are then assigned to each possible process and are 
drawn via Monte Carlo. Although process (1) produces a 
free electron and process (3) changes the energy of a free 
electron, we do not trace them in our simulations, but 
instead simply assume that they are thermalized before 
any further interaction takes place. 

As mentioned earlier, we assume that photoioniza- 
tion equilibrium is established locally everywhere in the 
plasma. Therefore, photoionization of H- or He-like ions 
is always followed by radiative recombination and radia- 
tive cascades to the ground level. For ions of lower charge 
states, K-shell photoionization is followed by either fluo- 
rescent emission or Auger decay. After photoexcitation, 
the ion in an excited state produces one or more pho- 
tons in the downward radiative transitions that eventu- 
ally lead to the ground level. By the combination of the 
photoexcitation and the downward radiative transitions, 
the resonance scattering of the lines is incorporated in 
the Monte Carlo simulation. 

For the H-like and He-like ions, we compute all of the 
atomic quantities with the Flexible Atomic Code (FAC) 9 . 
This includes energy levels, oscillator strengths, radia- 
tive decay rates, photoionization cross sections from all 
levels between the ground state up to n = 10. Radia- 
tive recombination rates onto each of levels as a function 
of the electron temperature are also computed from the 
photoionization cross sections using the Milne relation. 
Collisional transfers and UV photoexcitations in the He- 
like ions are not included in the current version of the 
code. 

Photoionization of ions with three or more electrons 
de-excites by either an emission of fluorescent photon or 
ejection of one or more Auger electrons. The K-shell pho- 
toionization cross sect ions are computed using the fitting 
formulae provided by iBand et alJ l)1990l) , which roughly 
accounts for the change in the K-edge energy and cross 
section as functions of the charge state. The fluorescence 
yields and Ka energies, however, are assumed to be fixed 
for a given ele ment and their values are adopted from 
(iLarkmsl 119771) and l|Salem. Panossian fe Krausd Il974t 
!Krausell979l) . respectively, for neutral atoms. We also do 
not account for K-shell ionization of Li-like ions, which 



8 http://hcasarc.gsfc.nasa.gov/docs/softwarc/xstar/xstar.html 9 http://kipac-tree.stanford.edu/fac/ 
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can lead to the production of the He-like forbidden line. 
We use L-shell cross sections from EPDL97 10 , which is 
also valid for neutral atoms. The Li, L2 and L3-shell 
cross sections are taken into account for the ions with 
three, five, and six electrons, respectively. Finally, we do 
not include L-shell fluorescent emission, since the yields 
are typically much smaller compared to the K-shell fluo- 
rescence yields. 

We consider Compton scattering only by free electrons, 
which is an excellent approximation for highly ionized 
media. The differential cross section is given by the 
the Klein-Nishina formula and the electrons velocities 
are drawn from a Maxwellian distribution with the lo- 
cal electron temperature. 

Doppler shifts due to the motion of the stellar wind are 
accounted for in all of the processes. When a photon en- 
ters a cell, the cross sections are calculated for the photon 
energy in the co-moving frame and Lorentz transformed 
back into the rest frame. The Doppler shifts could affect 
the line strength through the resonance line scattering. A 
velocity gradient can change the ratio of the line escape. 
In this Monte Carlo simulation, such Doppler effects on 
the lines are included, though the reproducibility of the 
velocity gradient is limited by the geometry cell size. 

4. THE PHYSICAL PROPERTIES OF THE STELLAR WIND 
IN VELA X-l 

4.1. Structure of the Stellar Wind 

We first search for the value of the mass loss rate that 
reproduces the observed line intensities as a function of 
the orbital phase. We assume that the velocity structure 
of the stellar wind is represent ed by a generalized CAK 
(ICastor. Abbott fc Kleinlll97^ profile: 

v(r) = Voo (l - ^ , (1) 

where r is the distance from the center of the star, i?„ is 
the stellar radius, and is the terminal velocity. We use 
a terminal velocity of 1100 km s" 1 as determined from 
UV observations of Vela X-l fP rinia etafl Il990tl and 
fix th e value of to 0.80 l)Pauldrach. Puis fc Kudritzkl 
EH). 

For a given velocity profile and a mass loss rate M* , the 
wind density can be calculated by applying the equation 
of mass continuity assuming spherical symmetry, 



■iit /im p v(r)r 2 ' 

where /i is the gas mass per hydrogen atom, which is 
~ 1.3 for cosmic chemical abundances. This specifies 
the density at each point in the stellar wind, and the 
observed X-ray luminosity can be used to obtain a map 
of the ionization parameter. 

The binary system parameters used in the simulations 
are listed in Table The X-ray radiation field from 
the neutron star is modeled based on the parameter de- 
termined from the fits to the observed spectra at phase 
0.25 and phase 0.5. Since no change in the average in- 
trinsic luminosity is observed between these phases, we 
assume that the luminosity remains unchanged during 

10 http://www.llnl.gov/cullenl/photon.htm, and distributed 
with Geant4 low energy electromagnetic package 



the eclipse phase as well. The power-law spectrum with 
a photon index r = 1 is assumed to extend up to 20 keV, 
which is the cut-off energy det ected by previous hard 
X-ray obser vations with Ginqa ( Makishi ma et al1ll999(l 
and RXTE ijKrevkenbohm et al.lll999j) . The correspond- 
ing 0.5-20 keV luminosity is 3.5 x 10 36 erg s _1 . In the 
following subsection, we discuss the results of our sim- 
ulation for mass loss rates of 5.0 x 10~ 7 , 1.0 x 10~ 6 , 
1.5 x 10~ 6 , and 2.0 x 10" 6 M yr -1 . 

4.2. The Ionization Structure 

Figure ITT1 shows maps of the H-like ion fraction of Si 
/siH-iikc at three different values for the mass loss rate. 
At a mass loss rate of 5.0 x 10~ 7 M© yr -1 , the concen- 
tration of H-like Si is highest near the mid-plane of the 
binary system as shown on the left panel of Fig. ^2 As 
the mass loss rate is increased, the density everywhere in 
the wind is also increased. For a fixed ionization param- 
eter, the distance to the neutron star has to be reduced, 
so the region moves closer towards the neutron star as 
shown in the middle and right panels of Fig. ^2 Mul- 
tiplying the maps in Fig. ^] by the local density of the 
stellar winds gives the maps in Fig. H3 These show the 
absolute abundances of H-like Si nsiH-Hke and the re- 
gion where the X-ray photons of the emission line are 
produced. 

According to these simulations, a large fraction of emis- 
sion lines from highly ionized ions such as H-like Si are 
produced in the region between the companion star and 
the neutron star. This is the region that is obscured by 
the companion star during eclipse. It turns out that the 
mass loss rate is sensitive to the ratio of line fluxes be- 
tween phase 0.5 and eclipse, because the size of the region 
becomes larger when we increase the mass loss rate. If a 
different value for the terminal velocity is assumed, the 
mass loss rate will differ such that the quantity M^/voo 
remains fixed. 

4.3. Estimate of the Mass Loss Rate 

The simulations are run for the various values of the 
mass loss rate. In order to calculate the observed spectra 
at different orbital phases, we designate the direction of 
the outgoing photon by the angle 9 that the momentum 
vector makes relative to the line of centers (see, Fig. 11311 . 
The phase 0.5 spectrum is constructed from the outgo- 
ing photons with 9 < 10° and the eclipse spectrum is 
made from photons with 170° < 9 < 180°. In this cal- 
culation, photons directly coming from the neutron star 
are neglected since only the emission lines and scattered 
components are of interest. 

Although many emission lines are resolved in the Chan- 
dra spectra, we use only the Si XIV Lya line to estimate 
the mass loss rate. This line is one of the brightest lines 
observed in the data. Since the energy of this line is suffi- 
ciently high, the uncertainties associated with foreground 
absorption is lower than the low energy lines. At this en- 
ergy, an uncertainty of 1 x 10 21 cm -2 in iVn, corresponds 
to only a ~ 5% uncertainty in the Lya intensity. 

The simulated intensity ratio of the Si XIV Lya lin 
between phase 0.5 and eclipse as a function of the mass 
loss rate is shown in Fig. [21 The errors associated with 
these points are due to the finite number of photons gen- 
crated in the simulation. The horizontal lines indicate 
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the lex range in the observed value. As mentioned earlier 
in § 14.21 the ratio becomes larger as the mass loss rate 
is increased. The intensities of the Si Lya in the phase 
0.5 and in the eclipse are shown in Fig. El From these 
three plots in Fig. El an d Fig. El the mass loss rate is 
estimated to be (1.5-2.0) xlO -6 M©yr -1 . 

The estimated mass loss rate of the stellar wind is con- 
sistent with the observed X-ray continuum luminosity. 
The mass accreting rate onto the neutron star M acc is 
given by 



lUacc - 4E)2 , 



(3) 



where M* is the mass loss rate of the stellar wind, and 
D is the distance of the neutron star from the center of 
the OB star. R acc is the accretion radius express as: 

2GM ns 



Ra 



(4) 



where v re i is the relative velocity of the neutron star 
and the stellar wind. The gravitational energy of the 
accreting material is converted into X-rays. The X-ray 
luminosity resulting from this accretion will simply be 
the rate at which gravitational energy is released: 



L, 



GAL 



(GM ns ) M* 



R 



(5) 



where it is assumed that most of this energy is liber- 
ated near the neutron star surface (of radius R nR ). If 



M nK = 1.7M Ijvan Paradiis et alJ Il977t iBarziv et all 
I2001|) . R ns = 10 km, v rel = 640 km s" 1 (v wmd = 
570 km s" 1 , Vorut = 300 km s" 1 ), D = 53.4i? Q and 
M* = 1.5 x 10 _6 M©yr _1 are applied, an X-ray luminos- 
ity of Lx = 4.7 x 10 36 erg s _1 is obtained, which agrees 
with the observed luminosity of 3.5 x 10 36 erg s _1 within 
a factor of 2. 

From the ecl ipse spectrum obtained by ASCA, 
iSako et ail l)1999j) estimated that the mass loss rate as- 
sociated with the highly ionized component of the stellar 
wind is 2.7 x 10 _7 M Q yr~ 1 , which is a factor of 5-8 lower 
than the current estimate. The lower mass loss rate was 
derived most likely because only the eclipse data was 
used. As the mass loss rate is increased, the intensities 
of the emission lines from highly ionized ions become 
higher. When the mass loss rate is increased further, the 
photoionized region is occulted by the companion star 
and the intensities of the emission lines during the eclipse 
phase become lower. This tendency for the Si XIV is 
shown in Fig. 1151 (left) . There were two solutions for the 
eclipse spectrum, and, ISako et all ilHHHti found the solu- 
tion of the lower mass rate. However, we note that the 
intensity of emission lines in the phase 0.50 obtained with 
Chandra cannot be explained with such a lower mass 
loss rate. We, therefore, have concluded the higher mass 
loss rate of 1.5-2.0 xlO -6 M0yr -1 . These measurements 
were probably not possible with the moderate spectral 
resolution of the instruments on board ASCA. 

4.4. The Global X-ray Spectrum 

Adopting a mass loss rate of 1.5 x 10~ 6 Mq yr _1 es- 
timated from the Si XIV line, we use the results of the 
Monte Carlo simulation to generate broadband spectra 
including both the continuum and emission lines from 



other ionic species for the three orbital phases observed 
with the HETGS. Photons with outgoing momenta in the 
6 ranges of 0°-10°, 85°-95°, and 170°-180° are used for 
the phase 0.50, phase 0.25, and eclipse model spectra, 
respectively. 

Figure El shows both the observed spectra above 
~ 2 keV and the simulated model spectra convolved 
through the response of the HETGS. The simulated spec- 
tra are multiplied by the effects of inter-stellar absorp- 
tion (Ah = 6 x 10 21 cm~ 2 ) in addition to intrinsic ab- 
sorption by the stellar wind, which depends on the ob- 
served binary phase. This matches the observed con- 
tinuum shape remarkably well at both phases 0.25 and 
eclipse. At phase 0.5, however, the observed column den- 
sity of Ah = 1.85 x 10 23 cm -2 is in excess of the Galactic 
plus stellar wind by an amount Ah = 1.7 x 10 23 cm~ 2 . 
We apply this additional absorption component to the 
continuum emission to match the data. 

The simulated spectra at the three orbital phases are 
shown in Fig. 1171 The observed emission line spectra in 
comparison with the simulations are shown in Fig. 1181 
The lack of obvious discrete spectral features at phase 
0.25 is consistent with the model predictions, which show 
P-Cygni lines with comparable emission and absorption 
equivalent widths and are therefore washed out after con- 
volution with the spectral response. We note that the 
additional absorption component observed at phase 0.50 
does not appear to be affecting the emission lines, which 
implies that this absorber is localized to a small region 
along the line of centers and in front of the neutron star 
at phase 0.5. The simulated spectra match fairly well 
the emission lines from H-like and He-like ions, as well 
as the radiative recombination continuum emission. The 
ratios of simulated line intensities to the observed ones 
are plotted for those lines that are detected at greater 
than 5 a in Fig. El All of the simulated line inten- 
sities are within the observed values to within a factor 
of ~ 3. One of the uncertainties of the line intensities 
is the absorption, and it is difficult to derive the line 
intensities precisely by unfolding the absorption effect. 
In the simulation, the effects of inter-stellar absorption 
(Ah = 6 x 10 21 cm -2 ) and intrinsic absorption by the 
stellar wind (Ah ~ 1.5 x 10 22 cm' 2 ) are assumed. If 
there is an additional absorption of 10% of the assumed 
value, the intensities of Ne lines are reduced to ~ 50%. 
Additionally, the absorption cross sections are very sen- 
sitive to the ionization degree of the absorption matter. 

4.5. The Iron Ka Line Complex 

Fluorescence emission lines, in particular the bright 
iron line complex, provide information about the distri- 
bution of cold material around the X-ray source. Fig- 
ure |S1 shows the iron Ka line spectrum at the three or- 
bital phases, together with the models obtained from the 
simulations described in the previous subsection. During 
eclipse, the iron Ka line intensity calculated from the 
simulation is consistent with the observed one. There- 
fore, it seems likely that the iron Ka line observed in 
eclipse is produced in the cool regions of the extended 
stellar wind. On the other hand, we find that this same 
model cannot explain the observed intensity of the Ka 
line in the other two phases. The equivalent widths from 
the simulations are only 4 eV and 11 eV for phases 0.25 
and 0.50, whereas the observed values are 51 eV and 
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116 eV, respectively. 

As an obvious source of additional line production re- 
gions, we first consider the surface of the companion star. 
Using Monte Carlo simulations and adopting the geomet- 
rical parameters listed in Table we obtain reflected 
iron Ka equivalent widths of 29 eV and 65 eV during 
phase 0.25 and phase 0.50, respectively, assuming that 
the metal abundance of the companion star is also 0.75 
times the cosmic value. Therefore, the stellar wind and 
the companion surface combined yield equivalent widths 
of 33 eV (phase 0.25) and 76 eV (phase 0.50). The ob- 
served values are still in excess these values by 18 eV and 
40 eV at phase 0.25 and 0.5, respectively. 

We note the presence of excess absorption of 
Nh = 1.7 x 10 23 cm -2 observed during phase 0.5 rel- 
ative to that at phase 0.25. This suggests that a distri- 
bution of cold localized material lies along the line of cen- 
ters somewhere behind the neutron star as viewed from 
the companion. We therefore expect that at least some 
amount of Fe-K line photons are pro duced in this cloud. 
This has already been pointed out bv llnoua l)1985fl using 
the Tenma data. 

In order to estimate the contribution of the cloud to 
the iron line flux, we have performed a Monte Carlo simu- 
lation of a partially covered cloud irradiated by an X-ray 
source. The geometrical configuration is shown schemat- 
ically in Fig. [201 The thickness of the cloud is assumed 
to be TVjj = 1.7 x 10 23 cm~ 2 and we assume a uniform 
density. The cloud is assumed to be neutral and its metal 
abundance is set to 0.75 cosmic. The resulting emission 
spectra as functions of the orbital phase with the solid 
angle as the free parameter are computed through the 
same procedures as in the previous simulations, and the 
simulated Fe Ka line equivalent widths are measured at 
the three orbital phases. 

Figure |^ shows the relation between the solid an- 
gle subtended by the cloud and the iron line equiva- 
lent width. As can be seen from this figure, the missing 
equivalent widths outside of eclipse (18 eV at phase 0.25 
and 40 eV at phase 0.50) can be explained if the cloud 
covers 25-40 % of the sky as viewed from the neutron 
star. Although this geometrical configuration may be 
oversimplified, its general characteristics are consistent 
with the fact that the spectrum observed at phase 0.25 
does not exhibit heavy absorption and the fa ct that soft 
X-ray lines are still observed during phase 0.5. iSato et aD 
( 1986) reported a gradual increase of average absorption 
column density between orbital phase ~ 0.2 and ~ 0.9 
from iV H ~ 1 x 10 22 cm" 2 to N H ~ 3 x 10 23 cm" 2 
observed in the Tenma data of Vela X-l. 

The enhanced absorption at orbital phase 0.50, and 
the localized cloud responsible for the iron line produc- 
tion can be attributed to an "accretion wake" , which is 
produced by a combination of several physical effects, 
including gravitational, rotational, and radiation pres- 
sur e forces, and X-ray heating. Numer ic al simulations 
by Sawada. Matsuda & Hachisu (1986), Blondin et al. 
( 1990), Blondi m Stevens fc Kallmanl l|1991fl and Blondin 
(1994) exhibit such structures trailing the accreting neu- 
tron star. 

4.6. The Velocity Field of the Stellar Wind 
4.6.1. Mismatch between observation and simulation 



When the observed line profiles of a few of the brightest 
lines are compared with the output of the Monte Carlo 
simulations, wc find that there is a mismatch in the aver- 
age velocity shift at both phase 0.5 and eclipse. FigureEl 
shows comparison of the observed Si XIV and Mg XII 
Lya lines with the simulations convolved with the in- 
strument response. Although the sign of the shifts are 
correct, the simulations seem to predict a larger velocity 
shift than seen in in the data. Figure[3shows the velocity 
shifts for all of the bright highly ionized emission lines 
detected in the spectrum. The average velocity shifts 
(Aw) observed between phase 0.5 and eclipse lie range 
in ~ 300-600 km s _1 whereas the simulations predict 
Av ~ 1000 km s _1 , which is larger approximately by a 
factor of two. 

4.6.2. Interaction between X-rays and the stellar wind 

One of the most obvious candidate mechanisms that is 
capable of reducing the wind velocity is X-ray photoion- 
ization. In the CAK-model of isolated OB stars, a stellar 
wind is driven by the radiation force on the large num- 
ber of resonance lines that lie in the UV. In particular, 
the dominant ions that are responsible for this are the L- 
shell ions of C, N and O, which have strong valence shell 
transitions in the UV range and have large abundances 
compared to the other metals. 

In HMXBs, there are regions in the stellar wind where 
X-ray photoionization destroys the L-shell ions of C, N, 
and O into K-shell ions and, therefore, reduce the amount 
of radiative force. This effect is maximized near the neu- 
tron star where the X-ray flux is high and which also hap- 
pens to coincide with the region where the wind would 
otherwise be accelerated in the absence of the X-ray 
source. As a result, the stellar wind velocity within this 
region should be lower than that predicted by the CAK- 
model. The effect of X-ray irradiation on line-driven stel- 
lar win ds in HMXB systems has be en st udied theoreti - 
cally bv lMacCregor fc Vitellol Jl982) and lMasail Ip84^ . 
and they have also been seen in UV spectral observations 
ijKaner et alJlTflfl^ . 

4.6.3. One dimensional calculation of the velocity 
structure 

In order to investigate how the stellar wind velocity 
is affected due to X-ray ionization, we have performed 
a simple one-dimensional calculation to estimate the ve- 
locity profile along the line of centers. The extension of 
this calculation to three-dimensions is required to self- 
consistently incorporate into the Monte Carlo simulator 
and predict the emission line spectra, which is beyond 
the scope of this paper. Nevertheless, as we show be- 
low, this estimate in 1-D is still useful for understanding 
semi-quantitatively the mismatch between the observed 
and predicted velocity shifts. 

In this calculation, we make two simplifying assump- 
tions. We first assume that the wind is driven radiatively 
by UV resonance absorption and that there are no other 
relevant forces that provide additional acceleration. The 
other assumption is that the force is proportional to the 
local number densities of C, N and O ions with more than 
one electron in the L-shell (i.e., Li-like and lower charge 
states). Figure [23 shows the velocity as a function of dis- 
tance along the line of centers in the CAK-model and the 
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radiative force, which is numerically derived from the 
CAK-model velocity structure. As can been seen from 
these figures, the radiation force peaks in the region be- 
tween the companion star surface and the location of the 
neutron star. This is, however, the region where most of 
the X-ray lines are produced in this system. Therefore, 
the acceleration of the stellar wind in this region must 
be reduced since the C, N, and O ions inside this region 
are primarily H-likc, He-like, or fully-stripped, and they 
alone cannot provide the force required to maintain the 
acceleration. 

We estimate the magnitude of this effect as follows. We 
start with the CAKvelocity profile, which has an implicit 
assumption about the force as a function of distance, and 
can calculate the density everywhere in the wind. This 
can then be used to derive the ionization structure and, 
in particular, the distribution of C, N, and O ions, which 
provides an estimate of the new force and hence the ve- 
locity profile. This process can be repeated and iterated 
to find an ionization structure that is self-consistent with 
the velocity profile through the following steps, 

1. start from the velocity structure (v[(r) = Vo(r) 
(CAK-model)), and, calculate a force structure 
(/i(r)), then, obtain v\{r). 

2. start from v' 2 (r) = Vi(r), and, calculate /2(f); then 
obtain V2(r). 

3. start fromi>3(r) = (vi(r)+V2(r))/2, and, calculate 
/3(r), then obtain i»3(r). 

4. start from v' 4 (r) = t>3(r), and, calculate fi{r), then 
obtain v^r). 

5. start from 1)5(7-) = (^(r) + Vi{r))/2, and, calculate 
/5(r), then obtain v 5 (r). 

6. start from v' 6 (r) — v§(r), and, calculate f$(r), then 
obtain Ve(r). 

The computed velocity structures («o( r )> v 2(r), 
«3(r), «4(r), «5(r) and ve(r)) are plotted in Fig. [231 As 
shown, the solution has converged after the 5th itera- 
tion. In this simple model, the stellar wind velocity at 
the point of the neutron star is only ~ 180 km s _1 com- 
pared to ~ 570 km s _1 predicted by the CAK-model. 
This reduction in velocity of ~ 390 km s _1 is remarkably 
similar to the observed velocity offset between the model 
and data observed during phase 0.5, where most of the 
emission is believed to originate from near the neutron 
star. 

Our one dimensional calculations confirm that the pho- 
toionization by the X-ray radiation can make the stellar 
wind near the neutron star slower by a factor of 2-3 com- 
pared to an undisturbed CAK wind. This results are in 
agreement with the observed Doppler shifts in the emis- 
sion lines. The X-ray radiation of the neutron star is 
disturbing the flow of the CAK-model stellar wind in the 
Vela X-l system, and then, the velocity of the stellar 
wind in the highly photoionized region becomes slower 
than that estimated from CAK-model. 

In order to understand the flow of stellar wind precisely 
and obtain the X-ray spectrum from the photoionized 
stellar wind, further considerations should be needed. 
A three dimensional calculation of the electromagnetic 



hydrodynamics would be required for a fully quantita- 
tive analysis. The lower velocity of the stellar wind due 
to the X-ray photoionization leads to a density increase 
and changes of emission line intensities. In that case, an 
adjustment of the mass loss rate could be needed to re- 
produce the observed spectra. Though the mass loss rate 
is proportional to the density and the velocity, how large 
the adjustment is required is not simple. The neutron 
star goes around the companion star, and, the photoion- 
ized region moves in association with the orbital motion 
of the neutron star. It takes about 0.5 days for the stel- 
lar wind to reach the position of the neutron star from 
the surface of the companion star assuming the CAK- 
model. It is not negligible in comparison to the orbital 
period of the neutron star (8.96 days), and, calculations 
including the time evolution due to the orbital motion 
are required. Additionally, gas pressure forces not only 
in the radial directions but along the circumferential di- 
rections are important since the stellar wind flow is not 
spherically symmetry any more. 

5. SUMMARY AND CONCLUSIONS 

Vela X-l was observed by the Chandra HETGS at three 
distinct orbital phase ranges. In order to understand 
the geometry and the dynamics of the X-ray irradiated 
stellar wind in this system, we have developed a Monte 
Carlo spectral simulator to predict the emission line spec- 
trum for a given set of stellar wind parameters. The 
model spectra were then compared to the observed spec- 
tra to construct a self-consistent picture that explains 
the global properties of the X-ray spectra. The qualita- 
tive aspects of this model and its various components are 
illustrated in Figl2H 

The X-ray emission line intensities of highly ion- 
ized H- and He-like ions are well-reproduced by a 
CAK-model of a stellar wind with a mass loss rate 
of (1.5-2.0) xlO _6 M Q yr _1 and a terminal velocity of 
1100 km s _1 irradiated by the X-ray continuum radia- 
tion from the neutron star. Through comparisons of the 
spectra observed at eclipse and at phase 0.50, we find 
that most of the X-ray recombination emission lines are 
produced in the region between the neutron star and the 
companion star (see Fig. I24f) where both the gas density 
and the local X-ray flux are high. This region is occulted 
by the companion star during the eclipse phase, which 
explains the relative strengths of the emission line in- 
tensities during phase 0.50 and eclipse. The same set of 
model parameters also nicely explains the lack of obvious 
discrete spectral features observed at phase 0.25. 

The properties of the intense iron fluorescent lines, 
which are observed throughout the binary orbit, are sen- 
sitive to the distribution of cold material relative to the 
location of the neutron star. Although the iron fluo- 
rescent line intensity observed during eclipse can be ex- 
plained entirely by fluorescence in the extended stellar 
wind, the line intensities during the other two phases are 
far brighter than predicted by the same model. We find 
that, in addition to the stellar wind and the X-ray irradi- 
ated surface of the companion star, a cold cloud partially 
covering the neutron star is required to explain the data. 
From the quantitative analysis, we find that a cloud with 
a thickness of A^h = 1.7 x 10 23 cm~ 2 and covers 25-40 % 
of the solid angle viewed from the neutron star, explains 
both the iron Ka line intensities as well as the excess 
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TABLE 1 

Summary of Vela X-l Observations 



Label OBSID Start Date Orbital Phase Exposure (sec) 



0.25 1928 2001-02-05 05:29:55 0.237 - 0.278 29570 

0.50 1927 2001-02-07 09:57:17 0.481 - 0.522 29430 

eclipse 1926 2001-02-11 21:20:17 0.980 - 0.093 83150 



TABLE 2 

Properties of continuum spectra derived from spectral fits. 



Orbital 


JV H a 


Photon Index 


Observed Flux 


Luminosity' 


X 2 / d.o.f. 


Phase 


(10 22 cm- 2 ) 




(erg cm~ 2 s^ 1 ) 


(erg s" 1 ) 




0.25 


1.45 ± 0.03 


1.01 ± 0.01 


3.0 X 10" 9 


1.6 X 10 36 


2872./ 3318 


0.50 


18.5 ± 0.3 


1.01 (fixed) 


1.5 x 10~ 9 


1.6 x 10 36 


1000./ 1051 



Note. — Fitting regions are 1.0-10.0 keV and 3.0-10.0 keV from HEG for 0.25 and 
0.50 orbital phases, respectively. The iron K-line region (6.3-6.5 keV) is excluded. Errors 
correspond to 90 % confidence levels. 

a The metal abundance is assumed to be 0.75 cosmic. 
b 0. 5-10.0 keV luminosity corrected for absorption. 



column density observed during phase 0.5. 

The one observed feature that this model fails to re- 
produce is the velocity shift of the highly-ionized recom- 
bination emission lines, which are observed to be smaller 
in magnitude by several hundred km s^ 1 compared to 
the predictions. Using a simple 1-D model, we find that 
destruction of UV-absorbing ions by X-ray photoioniza- 
tion can explain this reduction in velocity. However, the 
fully self-consistent 3-D model including the X-ray pho- 
toionization effect still remains an issue to be resolved. 
A full 3-D electromagnetic hydrodynamics calculation in- 



cluding time dependent effects would be required for the 
further study. 

In addition, there arc much more things to be done 
with this data set. These include study of the Comp- 
ton shoulder in the iron K line and the ratio of He-like 
triplet lines, which constrain the physical conditions in 
the emission region. There is also much more to be im- 
proved in the Monte Carlo code. For example, the code 
for ions from Li-like to near neutral is needed for a de- 
tailed study of the Si K line complex shown in Fig. |3J 
Such topics should be subjects of the next paper. 
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TABLE 3 

Derived Parameters of emission lines in the 0.5 orbital phase spectrum 



Energy 


Sigma 


Flux a 


Candidate 


Line Shift 


(Key ) 


leV) 


f 1 ft — 5 rilintnn pm ^ t; ^ 




(km s _1 ) 


o.oyuo_ 0009 


4+ 3 - 4 


8.8+3.1 


Ca II-XII Ka 




2.622±0.001 


4 0+ 1 - 5 

4.u_ 1 2 


17.71 3 ; 2 


S XVI Lya 




2.46197+8-22222 


0.6+or 

— U.D 


6.811;? 


S XV r 




r\ oi i i Q + 0. 00072 

z.0111Z_q QQQ42 


2+ 19 
u -^— 2 


10.31 22 


S IV-VIII Ka 




2.00634±0. 00016 


1 .2+°- 2 


23.81?;* 


Si XIV Lya 


+ 127+24 


1.86614±0. 00022 


1.4±0.3 


13.6+?? 


Si XIII r 


+183+25 


1 pcr-jv+0. 00055 

I.O.JO.) 1 _(| , 


o.8±5:° 


o 1+0.7 
° — 0.8 


Si XIII i 


+257 + |? 


i.84i36+2 : 2222? 


2.2±0.3 


14.1 + 1.2 


Si XIII f 


+311+ 42 


1.74447±0. 00026 


2.4±0.3 


14.61?; 2 


Si II-VI Ka 




1 7?qqS+°- 00035 


0.4±g"I 


3.4+0.6 


Al XIII Lya? 




1 . 59900~^~ i\ uiiaro 

— U.UUUott 


0.7±J : ° 


1-3+81 

— (J. 4 


Al XII r ? 




1.57976±0. 00056 


2.0+0.6 


3.9+0.7 


Mg XI 




i W oi +0.00053 
i.oazoi_ 00056 


6+ - 8 


1 a+0.6 
1 - 8 -0.5 


Fe XXIV ? 


+IO2I29 


1.47282±0. 00014 


1.2±0.2 


23.5+1.8 


Mg XII Lya 


1 qc; 97 q+0. 00017 
±.OOZI»_q 00()1(i 


1.0±0.2 


14.7+1.5 


Mg XI r 


+120+U 


1 q4S4fi+ 00021 
1.04CS4D_ o 00 22 


0.7±0.3 


7 ,+ 1.2 


Mg XI i 


+801S 


1.33213±0. 00023 


1.0±0.3 


q i +1 - 4 
y.i_ 1 2 


Mg XI f 


+230l 5 i 


1 qnsns+O' 00035 
1.3U8US_„ ()004 _, 


1 1+ - 4 
L - L -o.s 


5.ol?;J 


Fe XXI ? 


1 07700 + 0. 0()()2f< 
l.Zi (OO_ 000 27 


1.1±0.3 


8.41?;* 


Ne X Ly7 




1.21160±0. 00021 


i o+ - 2 

l - u -0.3 


13.3l?;l 


Ne X Ly/3 




1 i 97-jo+0. 00048 
i.iz(oz_ () ()0()4<) 


i c+0.5 
1 -°-0.i 


7.41?;? 


Ne IX 




1 n74-?9+ ' 00018 
l.U^JZ_ 00035 


o.il8:! 


5.81?;? 


Ne IX 




1 02242+ ' 00011 
±.UZZ4Z_ 00012 


0-618;? 


44 2 +5.5 


Ne X Lya 


+1821 32 


0.92246±0. 00027 


0.8±0.3 


28.91^ 


Ne IX r 


+149+12. 


n m COK+0. 00033 
0.916251,, 00032 


1 4+0-3 


49.1+11 


Ne IX i 


+475l?°« 


u.yuooo_ 00055 


1 «+ - 5 

1-O-0.4 


37.i±i;i 


Ne IX f 


+165i? 6 8 ? 



Note. — Errors correspond to 90 % confidence level. 

a Inter stellar gas absorption is corrected. The hydrogen column density of 
6 X 10 21 cm~ 2 is assumed, corresponding to the density of 1 H cm~ 3 and the 
distance of 1.9 kpc. 
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TABLE 4 

Derived Parameters of emission lines in the eclipse phase spectrum 



Center Energy 


Sigma 


Flux a 


Candidate 


Line Shift 


(keV) 


(eV) 


(10" 6 photon cm" 2 s" 1 ) 




(km s- 1 ) 



,+0.0040 
L -0.0043 

2.95657+o.onii 9 

2.61857t~ 
2.31035±0.0010 
2.00339±0. 



3.69431] 



1.86299+.°,- 



1.85271 
1.83924 
1.74247 



1.65752+.° 
1.57719±°, 
1.470681° 
1.35057±0 
1.34238+^ 
1.33122tg 
1.30640+.°, 
1.27592±°, 
1.23618±0 
1.20980t° 
1. 072411° 
1.02075±o 
0.92077+.°, 
0.91493±0. 
0.90406+.°, 



00034 
00031 
00032 
00086 
00017 
00045 
00040 

20 

OOO.'iO 
0010 

ooo.v, 

57 

00000 
00025 
00024 
00031 
00069 

oooos 

40 

00050 
00061 
00066 
00044 
00045 
0011 
00053 
00050 
00035 
00030 
00020 
00019 
00039 
00038 
00052 
00028 
00027 



8 fi+ 4 - 3 

a.b_ 1 2 


Q q+3.0 

y.3_ 2 7 


Pa TT VTT XC r\, 
Vjo, 11— All IV OL 


o o+ 3 - 4 


4 g+2.4 


A r VT-TX Krv 


n fi+ 2 - 

°- b -0.6 


1 Q Q + 3.0 

ld.»_ 3 


Q Yin T , r rt , 

o Avi LtyQ 


2.2±1.3 


fi +4.7 


Q T\T \7TTT T^"--, 

b IV— Vlll Ka 


4. I ±U.4 


no q + 2.7 
zo.z_ 2 5 


oi aiv ijyct 


1.6+.0.3 


no /i +2.8 

23.4_ 26 


c; VTTT 

oi Alll r 


n n+11 

°-°-o.o 

o q+0.5 


9 s + 1.3 
2 - 8 -l.l 


ol Alll l 


91 1+ 2 - 5 

— 2 3 


Si XIII f 


1.7+0.3 


10 o+2.'4 
19.6_ 2 2 


C : T T "\ 7T T/" 

bi 11— VI Ka 


al -o.i 


9 a+1.2 
2 - 8 Io.9 


T\ ,f _ V"T TJV. VVTTT 

Mg Al or rC AA111 


i q+0.8 
1 - J -0.6 


D-0-1.4 


IVlg Al 


i 9 +0.3 


9fi q+3.6 


Mp- XTT T vrv 

1VJ.H, Al L Jjyu; 


1 5+°' 4 
1 -°-0.3 


25 8+ 3 ' 9 
zo -°-3.5 


Ms- XI r 


1.2±0.8 


6 3+ 2 ' 4 
— 2.0 


Mg XI i 


9 O+0.5 


20.3+jJ 


Mg XI f 


1 q+0.7 


c c+2.5 


Fe XXI ? 


n a+°- 5 

u - M -0.9 


8.7+fl 


Ne X Ly7 


1 6+ 1 ' 2 


5.9+Jj 


Fc XX ? 


1 q+0.5 




Ne X Ly/3 


0.8±0.4 


in 1+6.2 
ia ' i -5.3 


Ne IX 


1.0±0.2 




Ne X Lya 


9+°- 4 


53t?J 


Ne IX r 


1 1+ 04 


47+?° 


Ne IX i 


0.8±0.3 




Ne IX f 



-314±51 
-323+J , 

-173± 2 3 8 9 
"34t 7 3 5 



-334+.51 
-373±69 

+25tiJ° 4 



— 308±|| 
-400+.1H 

-333l°* 



Note. 



Errors correspond to 90 % confidence level. 



a Inter stellar gas absorption is corrected. The hydrogen column density of 6 X 10 21 cm 
is assumed, corresponding to the density of 1 H cm -3 and the distance of 1.9 kpc. 



TABLE 5 

Comparison between lines of the phase 0.50 and of the 
eclipse. 



line Flux(0.50)/Flux(cclipsc) Av (km s" 1 ) 



Si XIV Lya 
Si XIII resonance 
Si XIII forbidden 
Mg XII Lya 
Mg XI resonance 
Mg XI forbidden 
Ne X Lya 
Ne IX resonance 
Ne IX forbidden 



10.3 ± 1.3 
5.81 ± 0.83 
6.68 ± 0.95 
8.94 ± 1.37 
5.70 ± 1.00 
4.50 ± 1.01 
5.18 ± 1.04 
5.27 ± 2.14 
4.13 ± 1.49 



441 ± 56 
506 ± 62 
345 ± 86 
436 ± 58 
493 ± 78 
205 ± 123 
490 ± 60 
549 ± 153 
498 ± 198 



TABLE 6 

Derived Parameters of the Fe Ka line. 



Orbital Phase 


Energy 
(keV) 


Sigma 
(eV) 


Flux 

(10~ 4 photon cm" 2 1 


EW 
s-i) (eV) 


0.00 
0.25 
0.50 


6.3958±0.0022 

fi qqqo+0.0018 
u.o»»^_ ( , 0005 
fi qqfjr+O.OOll 
o-3aoo_ .0012 


7 9+3.9 

'• Z -5.3 
0+ 7 ' 4 

11 0+ 1 ' 8 
±J - U -1.9 


1.7±0.2 
19.2±\i 
34.011J 


844 
51 
116 



Note. - 



Errors correspond to 90 % confidence level. 
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TABLE 7 

Adopted parameters for the Vela X-l simulation 



Parameter 



Value Reference 



Geometry 

Binary Separation D 
Companion Star Radius i?» 
Stellar Wind • • • 
Velocity Structure v(r) 
Terminal Velocity tioo 

Metal Abundance (Z > 2) 

X-ray Radiation 
Luminosity(0.5-20 keV) 

Spectrum Shape 
Energy Range 



53.4 R Q van Kerkwiik et al. (1995) 
30.0 Rq van Kerkwiik et al. f!995) 



„(l-R*/rf 
1100 km s" 1 
0.80 
0.75 cosmic 



3.5 X 10 Jb erg s" 1 

Power-Law (r = 1.0) 
13.6 eV - 20.0 keV 



Castor. Abbott fc Klein (1975) 

Prinia et al. (1990) 

Pauldrach. Puis k, Kudritzki ( 19861 

lor typical OB-stars 

CBord et al. 19761 



This observation 
(1.6 X 10 36 erg s" 
This observation 



(0.5-10 keV)) 
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Fig. 1. — The orbit of the neutron star in the Vela X-l system. The observer is located along the horizontal axis on the right. The bold 
lines span the range of orbital phases covered by the Chandra observations. 
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Fig. 2. — The observed HEG spectra of Vela X-l. The green, blue, and red histograms show the spectra at phase 0.25, phase 0.50 and 
eclipse, respectively. The solid black lines represent best-fit continuum spectral models for phases 0.25 and 0.50 with parameters listed in 
Tabic 
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Fig. 3. — A detailed view of the Vela X-l spectrum at phase 0.50 - MEG data in red and HEG data in blue. The spectrum is dominated 
by discrete features below E ~ 3 keV. The vertical green lines label the observed emission lines that are listed in Tabic |3] 
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Fig. 4. — Same as in Figure[3]for the eclipse phase. Note the striking similarity of the global emission line spectrum compared to that 
seen at phase 0.50. The vertical green lines label the emission lines listed in Table HI 
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Fig. 5. — The spectra of the Si K line complex observed at each orbital phase - phase 0.25 (top), phase 0.50 (middle), eclipse (bottom). 
The phase 0.25 spectrum is dominated by the intense continuum with possible evidence for weak absorption lines. The sharp drop at 
E = 1.84 keV is due Si K photoelectric absorption in the detector. At phases 0.50 and eclipse, emission lines from a wide range of charge 
states are seen. 
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Fig. 6. — Line profiles of Si XIV Lyo in the HEG spectrum (left) and Mg XII Lya in the MEG spectrum (right) plotted at a function 
of the velocity. The blue and red histograms show the observed data at phase 0.50 and eclipse, respectively, and the smooth curves were 
made from the output of the Monte Carlo simulation (see § 4.6) convolved with the respective HEG/MEG response functions. Note the 
discrepancy between the data and model on the order of 200 — 400 km s _1 in the sense that the model predicts larger velocity amplitudes 
than the observed values. 
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Fig. 7. — The plot of the velocity shifts for emission lines from highly ionized ions. The circles represent values derived from the 
observations, and the squares are values derived from the Monte Carlo simulations (see § 4.6). Filled and open symbols are used to 
designate values derived from phase 0.50 and eclipse, respectively. As shown in Figurc|n]for two particular lines, the lines from the other 
ions also systematically show a smaller velocity amplitude compared to the simulations. 
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Fig. 8. — The blown-up spectra of the iron Ka region at the three orbital phases. The data points represent the HEG spectrum and the 
dark histograms show the Monte Carlo models, which include iron fluorescent X-rays from the stellar wind alone (see § 4.5). The model 
accounts for essentially all of the line flux observed during eclipse, but fails to reproduce the non-eclipse data by a large factor. 
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Fig. 9. — A schematic diagram showing the sequence of how the ionization structure is calculated in the Monte Carlo simulator. The 
continuum photons from the neutron star (left) is propagated through each spatial cell and transfered to the next cell accounting for 
absorption. 
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Fig. 10. — The Vela X-l geometry and spatial grids used in the calculation of the ionization structure (left) and the Monte Carlo spectral 
simulation (right). 
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Fig. 11. — The map of the H-like ion fraction of Si. 
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Fig. 12. — The density map of H-like Si. 
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Fig. 13. — The definition of the angle 9, the angle of the outgoing photon relative to the line of centers. 
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Fig. 14. — Intensity ratio of Si Lya between phase 0.5 and eclipse as a function of the mass-loss rate. Each point represents results 
derived from the Monte Carlo simulations, and the error bars corresponds to lcr Poisson error due to the limited number of trials. The 
two horizontal lines denote the lcr range of the observed value. This figure shows that the mass-loss rate in the Vela X-l system is 
~ 1.5 - 2.0 X 10~ 6 Mq yr- 1 . 
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Fig. 15. — The Si Lya line intensity during eclipse (left) and for pha se 0.50 (right) as a function of the mass-loss rate. The horizontal 
lines denote the la ranges of the observed values. As in Figure 1141 a consistent solution can be obtained for a mass-loss rate of ~ 
1.5 - 2.0 X 10" 6 M Q yr" 1 . 




Fig. 16. — Comparison of the simulated spectra with the observed data above 2 keV. The normalization of the simulation models are 
fixed between the three phases ranges, assuming that the average X-ray luminosity does not change. The simulated models are multiplied 
by a Galactic column density of Nu = 6 x 10 cm -2 and convolved with the instrument response. 
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Fig. 17. — Monte Carlo simulation spectra at the three phases for the mass loss rate of 1.5 X 10 Mq yr 1 . These spectra are not 
convolved with the instrumental response. 
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Fig. 18. — Emission line spectra obtained with Chandra MEG. The red curves represent the data and the green lines show the Monte 
Carlo data convolved with MEG instrument response. 
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Fig. 19. — The plots of the intensity ratios of the simulated lines and the observed lines. The filled 
from phase 0.50 and eclipse, respectively. Each of the measured lines intensities are within a factor of 3 



and open circles show the results 
from the predicted values. 
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Fig. 20. — Geometry of the cold cloud used in the Monte Carlo simulation to calculate the iron Koj equivalent width. 
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Fig. 21. — Relation between the cloud's solid angle and the equivalent width of the iron Ka line calculated with a Monte Carlo simulation. 
The thickness of the cloud is Ah = 1.7 X 10 23 cm -2 and assumed to be at a uniform density. Open and filled circles show the predicted 
equivalent width as observed at phases 0.25 and 0.50, respectively. 
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Fig. 22. — The velocity profile of the CAK-model stellar wind with Voc = 1100 km s _1 and (3 = 0.8 (left) and the force profile (right). 
The vertical dotted lines designate the location of the stellar surface and the neutron star. The force is numerically derived from the 
CAK-model velocity structure. Note that the force is most dominant in the region between the neutron star and the companion stellar 
surface. 
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Fig. 23. — The results of a 1-D calculation of the ste llar w ind velocity along the axis between the companion star and the neutron star 
accounting for X-ray photoionization as described in § 14.6.31 The mass loss rate is assumed to be 1.5 X 10 — 6 MQyr — 1 . The flattening of 
the velocity profile is due to X-ray photoionization destroying the ions that are capable of absorbing in the UV. 
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Fig. 24. — A conceptual picture of the Vela X-l system obtained from the present analysis. 



